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A semi-analytic theory for the pulse dynamics in similariton amplifiers and lasers is presented, 
based on a model pulse with adaptive shape. By changing a single parameter, this test function can 
be continuously tweaked between a pure Gaussian and a pure parabolic profile, and can even repre- 
sent sech-like pulses, the shape of a soliton. This approach allows us to describe the pulse evolution 
in the self-similar and other regimes of optical propagation. Employing the method of moments, the 
evolution equations for the characteristic pulse parameters are derived from the governing nonlinear 
Schrodinger/Ginzburg-Landau equation. Due to its greatly reduced complexity, this description al- 
lows for extensive parameter optimization, and can aid intuitive understanding of the dynamics. As 
an application of this approach, we model a soliton-similariton laser and validate the results against 
numerical simulations. This constitutes a semi-analytic model of the soliton-similariton laser. Due 
to the versatility of the model pulse, it can also prove useful in other application areas. 



I. INTRODUCTION 

Self-similarity is a recurring theme in strongly nonlin- 
ear systems. Its observation can be particularly infor- 
mative as it implies an underlying symmetry, which can 
be exploited mathematically through symmetry reduc- 
tion techniques In nonlinear optics, self-similarity 
emerges in the formation of Cantor-set fractals in mate- 
rials that support spatial solitons '5], the self-collapse of 
beams at high povifcrs [3] , and in the propagation of ultra- 
fast pulses of light in optical fiber amplifiers in the pres- 
ence of strong Kerr nonlinearity [4, 5]. In recent years, it 
was reported that self-similar propagation of short pulses 
in laser resonators is possible [g, 0] . These pulses have a 
nearly parabolic intensity profile and evolve self-similarly 
within the nonlinear segments of the laser cavity. Fiber 
lasers supporting self-similarly evolving pulses is now rec- 
ognized as new regime of pulse formation in the cavity 
of an ultrafast laser. This method is differentiated from 
the well-known solitary stretched-pulse (dispersion- 
managed) iQi] and all-normal-dispersion [10,] solutions to 
the Haus Master equation [ll| . There are interesting sim- 
ilarities as well as important differences between these 
regimes. From a practical point of view, the demonstra- 
tion of the similariton laser has led to the development of 
fiber lasers with significantly higher pulse energies fia . 
These fiber lasers are being studied by many groups [13l - 
[TgI , motivated by the various applications ultrafast lasers 
have in diverse areas of physics, from optical frequency 
metrology and material processing to next-generation ac- 
celerators. More recently, a new mode-locking regime, 
the soliton-similariton laser was reported, in which the 
pulse evolution is in the form of periodic alteration be- 
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tween soliton and similariton evolution [17[. One aspect 
of this regime is that the evolution is strongly nonlinear at 
every point in the laser cavity. The possibilities and lim- 
itations of this regime are largely in need of exploration, 
for which theoretical modeling is crucial. For all of these 
reasons, there is much desire to understand the physics 
of amplifier similaritons and self-similar lasers better. 

Numerical simulations provide good agreement with 
experiments [gI. [T^, [Tt} . However, they are computation- 
ally expensive, rendering extended explorations of the 
parameter space impractical. Moreover, a theoretical 
description can aid intuitive understanding of the dy- 
namics of self-similar evolution in optical amplifiers and 
lasers. Exact self-similar solutions have been derived for 
the optical pulse propagation in fibers with and with- 
out gain d, 1^, However, the pulse shape evolves 
during propagation, and the self-similar parabolic pulse 
profile is only asymptotically reached. Thus, several ap- 
proaches have been explored to derive a simplified de- 
scription which still captures the rich pulse dynamics in 
such systems. Based on various analytical methods, the 
pulse formation, pulse stability and energy scalability of 
similariton and other high-energy fiber lasers has been 
studied [19-;21.J. Also semi-analytic approaches, widely 
used in optics to investigate pulse propagation, have been 
employed. They aim to extract evolution equations for 
characteristic pulse parameters, reducing the partial dif- 
ferential equation for pulse propagation to a coupled set 
of ordinary differential equations. Such approaches are 
typically based on the method of moments (MOM) or a 
variational formalism, which have both been used to in- 
vestigate the evolution of the pulse energy and the tem- 
poral and spectral pulse width in the strongly nonlin- 
ear regime |22h24| |. Such studies typically rely on fixed 
pulse shapes such as Gaussian or sech pulses, yielding 
reasonable estimates for the pulse energy and duration, 
but no pulse shape information at all. An exception can 
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be found in '25|, where an adaptive super-Gaussian test 
function was used to investigate changes of the pulse pro- 
file during propagation. 

Here, we report on a semi-analytic theory for the pulse 
dynamics in similariton amplifiers and lasers including 
the soliton-similariton laser, based on a novel model pulse 
with adaptive shape. The key in this formulation is our 
ansatz function that can describe any pulse shape from 
a pure Gaussian to a pure parabolic profile, even includ- 
ing sech-like pulses (i.e., with sech^ intensity profile), the 
shape of a soliton. The pulse profile is tweaked by a sin- 
gle parameter, which is complemented by an additional 
degree of freedom for the pulse phase. This allows us to 
represent various pulse profiles as well as complex spec- 
tral shapes. Thus, our theoretical treatment appears to 
be capable of describing not only the self-similar but the 
other regimes as well, opening the way to a simple unified 
theoretical approach. 

Employing the method of moments [1^ [131 , the par- 
tial differential equation governing the pulse propagation 
is reduced to a finite set of coupled ordinary differential 
equations, which are much easier to analyze. In addition, 
the coefficients of the equations are helpful in forming an 
intuitive understanding of the dynamics by exposing the 
relative importance of the various effects. Through inves- 
tigation of these equations one gains access to valuable 
information about the pulse dynamics, e.g., of how ex- 
actly the various effects on the pulse are paired to balance 
each other to satisfy the periodic boundary conditions 
imposed by the laser resonator. Such information is ex- 
tremely difficult, if not impossible, to obtain by repeated 
numerical solutions of the full governing equation. Our 
approach is validated against numerical results for single- 
pass propagation and for the steady state dynamics of a 
soliton-similariton laser. 



II. TEST PULSE AND EVOLUTION 
EQUATIONS 

For propagation through a dispersive Kerr medium 
with a parabolic gain and instantaneously saturable ab- 
sorption, the evolution of the pulse envelope u{z,t) is de- 
scribed by a generalized nonlinear Schrodingerfor com- 
plex Ginzburg-Landau) equation of the form [28| 

\dzU — Dd^u + 7 u — \(^g + g,^d^ + r u. (1) 

Here, z and t are the propagation coordinate and the 
retarded time, respectively. D is the second order dis- 
persion coefficient, and 7 is the cubic nonlinearity pa- 
rameter. The dissipative processes are characterized by 
the central gain value g and spectral gain parameter g^ 
as well as the saturable absorption coefficient r. Gener- 
ally, D, 7, g, g^, and r are z dependent, since an optical 
system such as a fiber laser consists of a sequence of dif- 
ferent segments. Additionally, the parameter values can 
vary even within a segment, for example g if gain satu- 
ration is considered. 



A. Test Pulse 

For linear systems, 7 = r = 0, the complex Gaussian 

u (z, r) = A^/pi{t) exp {ipr'^ + 10) (2) 

with pi{t) ~ exp (— r'^) is an exact solution to Eq. ([T]), 
where r — t/T denotes the normalized time, and T{z), 
A(z), (t>{z) and l3{z) are the pulse duration, amplitude, 
phase and linear chirp parameter, respectively. Thus, 
for moderate nonlinearity, the Gaussian ansatz is still 
a good description of the steady-state pulse shape in a 
laser cavity 129143 In contrast, in the strongly non- 
linear limit, the pulse is approximately described by a 
self-similar pulse with a parabolic intensity profile. How- 
ever, an exactly parabolic pulse is an idealization and 
in practice the pulse shape is parabolic around the cen- 
ter, where most of the energy resides, but with a super- 
Gaussian fall-off in the wings [3| • Naturally, in the 
intermediate regime, the pulse shape combines features of 
a Gaussian pulse and a self-similar pulse. To reflect these 
properties, we have previously introduced a function of 
the type 

Pn (r) = exp r^Vfc j = 1 - + O {r'"+') (3) 

to describe the pulse profile, which represents a Gaussian 
for n = 1 and a parabolic profile for n ^ cxd [s^I- Here, 
the pulse duration T represents the Gaussian pulse width 
for n = 1 and half the total pulse width of a similariton 
for 71 -H- 00. This ansatz has been shown to be useful for 
the description of similariton lasers and trapped Bose- 
Einstein condensates 

A disadvantage of Eq. Q is that the pulse shape can- 
not be adapted continuously, but only in discrete steps. 
Using the Gauss hypergeometric function 2F1 for which 
efficient numerical evaluation routines exist (3^ , Eq. ([3]) 
can be expressed in closed form as 

Pn (r) = (l-r2)exp|^ [2F1 (l, n; 1 + n; r^) - l] | , 

(4) 

see also Appendix lAl In Eq. (jlj, n is not restricted to in- 
tegers, providing much more flexibility for describing dif- 
ferent pulse shapes. For example, sech^-like intensity pro- 
files, corresponding to a fundamental optical soliton, are 
very well represented by n « 0.5. Moreover, rather than 
a priory fixing n to a certain value, we allow n = n (z) to 
evolve during pulse propagation, describing the position 
dependent intensity profile together with the parameters 
A{z) and T {z). Along with n (z), the third order chirp 
parameter a (z) is introduced as a further degree of free- 
dom for the pulse phase in addition to (3 {z) and 0(z), 
to avoid mathematical problems with the evolution equa- 
tions for the pulse parameters [35|. The resulting ansatz 
for the envelope is given by 

u (z, r) = A^Jpn{T) exp (i/3r^ + iar^ + \(f)) . (5) 
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Naturally, Eq. Q is not the only function which is able to 
interpolate continuously between a parabolic and a Gaus- 
sian shape. In particular, the so-called q-Gaussian func- 
tion "SG*] has been used in various contexts, e.g., for the 
description of trapped Bose-Einstein condensates [37| . 
While the q-Gaussian has a somewhat simpler analyti- 
cal form, it is non-zero only on a finite interval (except 
for the limiting case of a Gaussian), which is unphysical 
for the applications considered in this paper. Addition- 
ally, our ansatz has the distinct advantage that it can also 
represent a sech^ profile to a very good approximation, 
which is essential for a versatile description of nonlinear 
optical propagation. 



B. Evolution Equations 

The generalized nonlinear Schrodinger equation Eq. ([1]) 
can be approximately solved by extracting evolution 
equations for the parameters of the model pulse Eq. ([S]). 
Here we use the method of moments (MoM) [2^; the 
derivation can be found in Appendix [BJ The resulting 
equations of motion are 
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where the prime denotes a partial derivative with respect 
to z. The weighing coefficients are given by 
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V(T)~M5rP„(T)]'dT. (11) 



To increase numerical efficiency, they are calculated only 
once for a sufficiently closely spaced n grid and tabulated. 
The evolution equations Eqs. ([6]) - (jlOp are also valid for 
z dependent coefficients in Eq. ([T]), which is especially 
important for effects like gain saturation. We note that 
the validity of the derived equations is not restricted to 
ansatz Eq. Q; in fact, they can be used for any such 
test pulse Pn with a continuously adjustable pulse shape 
parameter n (and p„ (r) = p„ (—■'')), like the q-Gaussian 
function [s^, Is^] ■ Only the weighing coefficients {n) , 
(n) and ijk (n) (Eq. (ITTt ') have then to be recalculated 
for that specific function. 



III. RESULTS 

To validate the ansatz Eq. ([S]), the equations of motion 
Eqs. - (jlOp are solved in different nonlinear propa- 
gation regimes. First, the soliton regime is considered. 
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characterized by negative dispersion and moderate non- 
linearity. Then, the self-similar propagation through gain 
fibers with positive dispersion is studied. Finally, the 
ansatz is employed to find the steady state solution of 
a soliton-similariton fiber laser, where alternate propa- 
gation in both regimes occurs. The equations of mo- 
tion Eqs. © - (jlOp are solved with a standard differen- 
tial equation solver, allowing for an efficient treatment 
of the problem. For comparison, also the results for 
the simplified Gaussian ansatz Eq. ([5^ are shown. The 
corresponding equations of motion [3l| can be obtained 
from Eqs. 0, © and ^ by setting n = 1, a = 
and n' = a' = 0. The semi-analytic results are vali- 
dated against exact analytical solutions of Eq. ([T]) or full 
numerical simulations, performed with a standard sym- 
metric split-step propagation algorithm [28i] . 



A. Fundamental Soliton 

For g = = r = 0, steady state solutions of 
Eq. ^ exist. For 7 > 0, £» < (or 7 < 0, Z? > 0), 
a special solution is given in form of the fundamental 



soliton, with the power 



A^sech^ (t/Ts), where 



Ts = A-^ (-2D/7)^/^ [281. To test the validity of our 
ansatz Eq. ([5]), wc extract the steady state solution of the 
evolution equations Eqs. (|6]) - ([T0| with g ^ g^j ~ r — 0, 
and compare it to the exact soliton solution. By setting 
dz = 0, we obtain /3 = a = 0, /i2'7o + SeoMo ~ 7?2/^o = 
which is fulfilled for n « 0.518, and tiol^^T^ = -?]oD. 
The pulse energy E — A^Teq can thus be written 
as S = eo{r1^/^i^^f'^A{~Dhf'^ « 2.7%A{-Dhf'\ 
The energy of the exact solution of Eq. ([1]), i.e., the funda- 
mental soliton, is — 2^/"^ A{—D /^)^^'^ , thus we have 
E « 0.99 i?s. The Gaussian ansatz, Eq. ([2]), is less ac- 
curate, yielding E « 1.05 i?s. In Fig.[l] the approximate 
(solid line) and exact (dashed line) solution are compared 
for a fixed pulse amplitude A. The results are virtually 
indistinguishable, demonstating that the ansatz Eq. ([S]) 
works very well in the soliton regime. For comparison, 
also the Gaussian steady state solution is displayed (dot- 
ted line). It provides a less accurate but still reasonable 
fit, even though it naturally fails to reproduce the char- 
acteristic sech^ soliton shape. 



B. Amplifier Similariton 

In order to test our ansatz in the self-similar regime, 
single-pass propagation in a gain fiber with positive dis- 
persion is studied. The investigated setup is as de- 
scribed in [4], with the fiber parameter values 7 — 
5.8 X IQ-^W-^mr^, D = 12.5 x lO^^ps^m^i, and 
g ~ 0.95 m^^; furthermore, r — g^ — Q. The initial pulse 
is assumed to be Gaussian {n — 1) with a fixed energy of 
12 pj. First, the pulse evolution is studied with ansatz 
Eq. ([SJ and by full numerical simulation for an initial 



pulse duration of 0.2 ps. Here, the pulse is characterized 
in terms of its temporal and spectral width TpwHM and 
/fwhm, respectively, which are the full width at half- 
maximum (FWHM) values of the instantaneous power 
and the power spectrum. Furthermore, n (z) is evaluated, 
describing the pulse shape of our ansatz Eq. ([S]) . For the 
numerical pulse, the kurtosis [13, Ull / {t — to)'^ p dt / a'^ 
is calculated, where p = P/ J Pdt is the normalized 
pulse power, to — Jtpdt = is the mean value, and 
= J (t — to)^ pdt is the variance; n it then extracted 
by determining the Pn in Eq. ^ with the same kurto- 
sis. In Fig. [21 the evolution of the pulse parameters is 
compared for the method of moments and full numeri- 
cal simulation. In Fig.[5Jc), s — n/ (n + l) rather than 
n itself is plotted to restrict the range of values to [0, 1]; 
i.e., s = 1/2 corresponds to a Gaussian and s = 1 to a 
parabolic pulse. In the example shown, s approaches 1, 
indicating that the pulse approaches self-similar evolu- 
tion. The overall agreement between semi-analytic and 
numerical results is excellent, indicating that our ap- 
proach works well also in the regime of self-similar propa- 
gation. Specifically, our ansatz Eq. ([5]) fully captures the 
transition of the pulse shape (see Fig.[5](c)). 

In Fig. 131 the instantaneous power and power spectrum 
are shown after a propagation distance of 3 m for Gaus- 
sian initial pulse widths (FWHM) of 0.1 ps (Fig.|3i;a), (b)) 
0.2 ps (Fig.lSIc), (d)) and 1 ps (Fig.E] (e), (f)), respec- 
tively. Ansatz Eq. ([5]) (solid lines) provides an excellent 
qualitative and quantitative approximation, reproducing 
very well the exact numerical pulse shapes and power 
spectra (dashed lines). The Gaussian approach (dotted 
lines) shows some deviations in the pulse duration and es- 
pecially the amplitude, but overall still provides a reason- 
able fit in time domain, see Fig.|3l[a), (c), (e). However, it 
naturally fails to reproduce the pulse shapes. Especially 
for strongly self-similar propagation as shown in Fig.[3Kc), 
where both our ansatz and the exact result exhibit a dis- 
tinct parabolic intensity profile, the Gaussian ansatz does 




Normalized time T/T 

s 

FIG. 1. Instantaneous power vs. time for the approximate 
and exact fundamental soliton solution; for comparison, also 
the Gaussian approximation is displayed. 
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not approximate the pulse shape well. Regarding the ob- 
tained power spectra, see Fig.[3][b), (d), (f), the Gaussian 
ansatz completely fails to reproduce the spectral features. 
The capability to faithfully reproduce spectral character- 
istics is particularly important from a practical point of 
view: Experimentally, optical spectra provide the most 
direct, immediately available and quite informative in- 
sight into the evolution of an ultrafast pulse among all 
the diagnostics at the disposal of the researcher. 



C. Soliton-Similariton Fiber Laser 

In the following, we apply our approach to self-similar 
propagation in a laser cavity, where the laser field is sub- 
ject to periodic boundary conditions in steady state op- 
eration. We choose a soliton-similariton laser setup as 
investigated in p7| , which is especially interesting in our 
context since the pulse undergoes self-similar propaga- 
tion as well as reshaping to Gaussian and sech^ profiles 
in the same cavity. In our case, the setup consists of a 
gain fiber, a piece of single mode fiber (SMF), a saturable 
absorber (SA), a bandpass filter, and again an SMF. The 
pulse evolves self-similarly in the gain fiber, and is tem- 
porally and spectrally filtered in the SA and bandpass 
filter, respectively. The group velocity dispersion (GVD) 
in the SMF is negative, approximately canceling the pos- 
itive GVD in the gain fiber. Several distinct nonlinear 
pulse shapes co-exist in the cavity: A parabolic profile is 
obtained towards the end of the the gain fiber, character- 
istic for self-similar evolution, then the pulse undergoes 
Gaussian spectral filtering and approaches a sech^ shape 
in the SMF, typical for a fundamental soliton. 

The parameter values for the gain fiber (SMF) are 
7 = 9.32 X lO-^W-^m-i (1.1 x lO-^W-^m-i), D = 
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FIG. 2. Evolution of the pulse duration, spectral width and 
pulse shape as a function of the propagation coordinate z, 
computed with the method of moments and by solving Eq. 
numerically. 



0.03845ps2m-i (-O.OlWps^m-i), on = 3.45m-i (0), 
and = 3.25 x 10^''ps^m~^ (0) [l3|. The gain is as- 
sumed to saturate with the pulse energy E, i.e., g = 
go/ {1 + E/Esnt), where i?sat — 2.21 nJ is the saturation 
energy. The bandpass filter is modeled by a segment of 
length L with g^L — 0.015 ps^, corresponding to a spec- 
tral width of 12 nm (FWHM), and the pulse power is 
additionally reduced by a factor of 5 to account for the 
overall linear loss of the optical cavity elements. For the 
SA, an unsaturated loss of go = 0-''' ^md a saturation 
power of Psat = 2.13 kW is assumed; its implementation 
is discussed in Appendix 

In Fig.m the MoM and full numerical results for the 
evolution of characteristic pulse parameters in the cavity 
are compared, where the sequence of optical elements and 
the fiber lengths are as indicated in Fig. Ufa). The pulse 
parameters are defined as described in Section llllBI The 
overall agreement between semi-analytic and numerical 
results is again excellent, compare Fig. Particularly, 
as can be seen in Fig.HJ^c), our ansatz Eq. ([S]) correctly 
predicts the almost parabolic pulse profile in the gain 
segment, with s = 1 for a parabolic pulse, the Gaus- 
sian shape after the filter (s = 1/2), as well as the sech^ 
shape in the SMF, corresponding to s « 1/3. In Fig.O 
the instantaneous power and power spectrum are shown 
after the gain fiber, before the SA, after the bandpass 
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FIG. 3. Instantaneous power and power spectrum, as ob- 
tained with the method of moments (solid lines), by full nu- 
merical simulations (dashed lines), and with the simplified 
Gaussian ansatz (dotted lines). The initial pulse durations 
are 0.1, 0.2 and Ips, respectively. 



6 



filter, and after the second SMF. The overah agreement 
between semi-analytic (solid lines) and numerical results 
(dashed lines) is very good both in the gain fiber and the 
SMF. Especially, ansatz Eq. ([S]) approximates well the 
distinct temporal and spectral pulse shapes in the differ- 
ent regimes. For comparison, also the Gaussian solution 
is displayed (dotted line). It provides a reasonable fit to 
the temporal and spectral width, but naturally cannot 
reproduce the pulse shape at all. Only after the band- 
pass filter, which forces the power spectrum to assume 
a Gaussian profile, the Gaussian ansatz closely matches 
the numerical solution (see Fig.[SKf), (h)). 



IV. CONCLUSION 



lent agreement. This constitutes a semi-analytic model 
for the soliton-similariton laser. A major advantage of 
the semi-analytic method is that the calculations are ap- 
proximately 100 times faster than the full numerical sim- 
ulations. This will allow the exploration of a vast pa- 
rameter range of interest to the design of fiber and solid 
state similariton lasers. Furthermore, this approach can 
be helpful for developing an intuitive understanding of 
the dynamics of self-similar evolution in optical fiber sys- 
tems by exposing the relative importance of the various 
effects. Due to the versatility of our test function, we 
expect it to also prove useful in other application areas 
in nonlinear optics, or in completely different fields such 
as the description of trapped Bose-Einstein condensates, 
as already exemplified in [33j . 



In conclusion, we have developed a semi-analytic the- 
ory for nonlinear optical ultrafast pulse propagation in 
the self-similar and other regimes, which we employ to 
study the pulse dynamics in similariton amplifiers and 
lasers. The key is the introduction of a model pulse 
with adaptive shape, which can continuously be tweaked 
with a single parameter to represent pulse shapes ranging 
from parabolic to Gaussian to sech^-like intensity pro- 
files. Thus, very different regimes of nonlinear optical 
propagation can be covered. Based on the method of 
moments, evolution equations are derived for the charac- 
teristic pulse parameters, specifying the pulse amplitude, 
duration, profile, and linear and third order chirp. Com- 
parison to exact analytical or full numerical results were 
performed for the soliton regime as well as similariton 
amplifiers and soliton-similariton lasers, showing excel- 
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FIG. 4. Evolution of the pulse duration, spectral width and 
pulse shape in the laser cavity, as obtained with the method 
of moments and by solving Eq. ^ numerically. 
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FIG. 5. Instantaneous power and power spectrum, as ob- 
tained with the method of moments (solid lines), by full nu- 
merical simulations (dashed lines), and with the simplified 
Gaussian ansatz (dotted lines). The results are shown at the 
positions i ((a), (b)), u ((c), (d)), iu ((e), (f)), and iv ((g), 
(h)) in the laser cavity, as indicated in Fig. Ha). 
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Appendix A: Test Pulse 

The test pulse Eq. ^ can be written as 
Pn (r) - (1 - t2) exp I ^ [2F1 (1, n; 1 + n; r^) - l] 



= (l-r2)exp{|r|^"[$ (r^ 1, n) - n-^] } , 

where 2F1 is the Gauss hypergeometric function and 
<i> is the Lerch Phi function, defined as $ (z, a, n ) — 
"^h^QZ^ I {n + k)°' for 1^1 < 1 and analytic contin- 
uation otherwise. For = 1, where 2F1 and $ 
both diverge, Pn has to be expressed in terms of 
the digamma function 4' (z) and Euler's constant 7, 
Pn (±1) = exp {--^ {n + l)- 7). 

These special functions are routinely implemented in 
many mathematical tools, and efhcient routines are avail- 
able [sj). However, we found it convenient to evaluate 
Eq. ([5]) by a series approach, using 



Pn (r) = (1 - T^) exp |t| 



2n 



1 2m 



m + n 



for < 1 and 



p„ (r) = (t2 - 1) exp |r 



2n \T\ 



-2m 



m>0 



cos(27rn) / n n 

X exp < TT — - — h TT 2 COS (Trn) — 1 tan —Tin 

sm(7rn) \2 

for > 1 (and n ^ N). For n e N, Pn is directly given 
by Eq. ©. 



Appendix B: Derivation of the Equations of Motion 

The equations of motion for the pulse parameters are 
derived using the method of moments [27)] . We intro- 
duce the energy Qq and the momentum Pg, 



Qo 



\u\ dt, 



-^0 = 2 / {u*tU ~ utu*)dt, 



Qi 



t\u\ dt, 



(t-io) |u| di, £> 1 



{t - toY (utu* - ulu) dt, £ > 



where to denotes the center of gravity. Due to the sym- 
metry properties of the ansatz Eq. ([5]), we have Qi = Q 
for odd i and = for even t, as well as to = 0. 

Multiplying Eq. ^ with u* and subtracting the com- 
plex conjugate, we can write 

\u\^ + Ddt (udtu* - u*dtu) = u*R - uR* , (Bl) 

with the dissipative term R — i (^g + g^df -\- r\uf'^ u. 

Multiplying with and integrating over t yields the equa- 
tions of motion for the Q^. Furthermore, multiplying 
Eq. ID) with ul and subtracting u* times the temporal 
derivative of Eq. ([1]) , and subsequently taking the real 
part of the resulting equation yields 



id, [ulu - utu*) - ADdt \ut\' + Ddf \u\^ - jdt \uf 



2 {utR* + u*R) - dt [uR* + u*R) . 



(B2) 



Multiplying with and integrating over t yields the equa- 
tions of motion for the Pi. We arrive at the evolution 
equations 



/oo 
{uR* - u*R) dt, 
-00 



(B3) 



/OO 
{uR* - u*R) dt, (B4) 
'OO 

/OO 
t'*(u*P-MP*)dt, (B5) 
-00 

/OO 
(-iD\ut\^ -jlufUt 
-OO' 

/OO nOO 
t (utR* + u;R) dt + i (uR* + u*R) dt, 
-00 J —00 

(B6) 

/OO 
\utf dt ~ 6DQ0 
-00 

/OO 
t^l^l^dt 
-00 

/OO pOO 
{utR* + ulR) dt + 3 / t^ {uR* + u*R) dt. 
-00 J —00 

(B7) 
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Inserting Eq. ([5]) , we obtain 

= 2geo + 2rA^^iQ 
1 



£0 2 



n'dneo 



9uT-^ ( - 2'/o - 8/^'^2 - 32a2£g _ 32/3a£4 ) , (B8) 



£2 (2— +3— j +n'a„e2 
= ~8DT~^ (/3£2 + 2ae4) + 2ge2 + 2rA^fi2 
+ 2g^T-^ ( -ir;2 + £o - 4/32^4 - IBa^es - 16/3a£6 



(B9) 



/A' T'\ 
£4(2-^ + 5—)+ n'a„£4 + 16DT-^ (/3£4 + 2a£6) 

= 25£4 + 2rylV4 + 9ojT-^[- \m + 12£2 - ^P^Eq 



32q; £10 — 32/3a£8 



(BIO) 



2^ + — j (/3£2 + 2a£4) + /3'£2 + /3n'9„£2 
2a'£4 + 2an' dnSA 



-DT-^ ( -?7o + 4/3^£2 + 16;3a£4 + 16a^£6 



7 ,2 



A>o + 25/3£2 + 45q;£4 + 2rA^ 13^12 + 4rA^Q;^4 
+ gc.r-2(^3/3£o - "^Pm + 42a£2 - 3ar/4 - ^%p'^ae& 
- 96;3a2£8 - 8/33£4 - 64a3£io) , (Ell) 



A' T' 

2— + 3— ) (/3£4 + 2a£6) + P'si + /3ri'9„£4 + 2a'£6 



2an'9„£6 + ^7AV2 + SDT" 



1 1 

2^0+4^2 



+ AjS'^EA + 16/3a£6 + iWes, 

= 2g/3£4 + 45a£6 + 2ryl^/3^4 + 4ryl^Q/i6 
3 



+ guT^"^ (21^62 - -Pm + 102a£4 - 3a?76 - 48/3^a£8 
- 96;3a2£io - 8/3^£6 - 64a3£i2y (B12) 



Eq. d?]) is obtained after multiplying Eq. (jB8|) by £2 /£o 
and subtracting Eq. (|B9p ; similarly, multiplying Eq. (jBSP 
by 3£2/£o and subtra cting Eq. (|B9t yields Eq. ([8]). Eq. dH) 
is obtained from Eq. (IB10[) by inserting Eqs. ([7]) and 
Furthermore, we derive Eq. (1101) by eliminating n'dnS2 
and n'dn£A from Eq. (|B11|) , using Eqs. (|B9|) and ([BT0| . 
respectively. Finally, Eq. © is derived from Eq. (|B12[) by 
eliminating (3n'dn£i with Eq. (jBlOl) and (3' with Eq. JTT 



Appendix C: Modeling of the Saturable Absorber 

In the Schrodinger equation Eq. ([ij , instantaneously 
saturable gain or loss is described by the term dzu\^^^ ~ 
r \u\^ u, with the solution 



u{L) 



Uo 



1 - 2rL\uo\ 



(CI) 



for an initial field uq and a propagation length L. Thus, 
the pulse power P (t) — \u {t)\ is transformed according 
to 



P(L) 



1 - 2rLPo ' 



(C2) 



while the phase of u is not altered. For saturable ab- 
sorption (r > 0), this approach only works in the weak 
field regime, i.e., 2rLPo <C 1. More generally, a saturable 
absorber (SA) can be modeled by the expression 



Pi = Poll 



1 + Po/Ps'. 



Po - qoP (L) , (C3) 



where qo is the unsaturated loss, and Pgat is the satura- 
tion power. 

In the following, we describe how to obtain the param- 
eter values of our test pulse Eq. ([5]) after passage through 
an SA of the form Eq. (|C3p . Most straightforwardly, 
this could be achieved by Taylor expansion of the pulse 
around its center at the in- and output of the SA and 
comparison of the leading terms \2§\. Here, we aim for 
a more global fitting method, consistent with the MoM. 
First, the equations of motion Eqs. ([S]) - (|5]) are solved 
for r = -1/ (2Psati) and g = g^ = D = j = 0, yielding 
the pulse parameters A (L), T (L) and n [L) of P [L) in 
Eq. (jC3p . The corresponding parameters Ai, Ti and ni 
for Pi are then derived by computing the 0th, 2nd and 
4th moment of Eq. ([C3t . 



t"'Pidt = / t"'Podt -qo t"'P [L) dt 

-00 J —00 J —oc 

(C4) 

with m = 0, 2 and 4, yielding 
i^m = A2r™+i£„. (no) - qoA^ (L) T"+i (L) Em (u (L)) 

- A2r™+i£„ (m) , (C5) 

with Sm defined in Eq. (jlip . From this, we obtain an 
implicit equation for ni, 



£0 (ni)£4 (ni) _ ypVi 
£2(ni) i/| 



and furthermore 



Ai = 



I So (ni) V2 
£2 (ni)i'o' 

TiEo {ni)' 



(C6) 

(C7) 
(C8) 
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The phase i/3 {t/Tf + ia [t/Tf + i0 of our test pulse ao {Ti/Taf. 
Eq. ([SJ is not altered, thus we get /3i = /3o (Ti/Tb)^, "i = 
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